MICROCOPY  RLSOlUl  ION  U bl  (MARI 

NA'K>NAl  HURfAU  ()»  ,!  ANIjAR!).,  i «.  ••  • 


THE  GENESIS  OF  DYNAMIC  SYSTEMS 
GOVERNED  BY  METZLER  MATRICES 


Technical  Report  No.  29 


Prepared  under  Contract  No.  00014-67-A-0298-0019 

Project  No.  NR-47-004 

for  the  Office  of  Naval  Research 


This  docuaent  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited. 

Reproduction  in  whole  or  part  is  permitted  for  any 
purpose  of  the  United  States  Government. 


Harvard  University 

1737  Cambridge  Street,  Room  404 

Cambridge,  Massachusetts  02138 


December,  1976 


THE  GENESIS  OF  DYNAMIC  SYSTEMS  GOVERNED  BY  METZLER  MATRICES 


0* 

Kenneth  J.  Arrow 

0.  Introduction. 

The  literature  on  dynamic  systems  In  economics  Is  vast,  and  an  Important  part 
of  that  deals  with  systems  of  differential  or  of  difference  equations  where  the 
Jacobian  of  the  right-hand  side  Is  a Metzler  matrix,  l.e. , a matrix  whose  off- 
diagonal  elements  are  non-negative.  Such  matrices  have  a wide  range  of 
applicability  In  dynamic  economic  models.  In  Input-output  analysis.  In  stability 
analysis  of  systems  of  excess  demands  governing  price  changes,  and  In  multi-sector 
and  multi-national  Keyt  >s1an  Income  determination  models.  Oskar  Morgenstern  early 
perceived  the  Importance  of  such  models  and  encouraged  research  In  them,  as  seen 
In  the  papers  by  Y.K.  Hong  and  M.A.  Woodbury  In  Morqenstem  (3).  For  a later 
survey,  see  (2). 

(The  bulk  of  this  work,  as  Indeed  the  bulk  of  the  work  on  dynamic  systems  In 
general,  concerns  what  might.  In  theological  terms,  be  called  the  eschatology  of 
' the  system,  the  questions  of  the  end  or  final  state  of  the  system.  In  this  paper, 

I want  to  concentrate  on  the  behavior  of  the  system  In  Its  Initial  phases.  Its 
genesis.  The  problems  revolve  mainly  about  the  presence  of  off-diagonal  zeroes 
In  the  matrix  governing  the  system  and  of  zeroes  In  some  components  of  the  Initial 
conditions  and  of  the  forcing  terms.  If  the  matrix  were  strictly  positive,  for 
example,  then  any  Initial  Impulse  anywhere  Immediately  (or  with  a lag  of  one  time 
period  In  a difference  equation  system)  produces  a positive  response  everywhere. 

But  If  there  are  some  zeroes  In  the  matrix,  then  the  transmission  of  the  Impulses 
Is  delayed.  In  the  case  of  differential  equation  systems,  the  effect  appears  as 
a lower  rate  of  growth  of  the  component;  Instead  of  Increasing  linearly.  It  may 
Increase  quadratlcally  or  even  with  a higher  power  from  Its  Initial  value  of  zero. 

The  discussion  of  this  subject  Involves  what  appears  to  be  a new  concept,  that 
of  first-positivity.  A sequence  Is  first-positive  If  the  first  non-zero  element 
is  positive.  Similar  definitions  can  be  given  for  sequences  of  vectors  and  matrices. 
The  methods  bring  together  a number  of  elementary  concepts  from  diverse  fields. 
Including  matrix  algebra,  the  theory  of  relations,  and  the  theory  of  differential 
equations. 

In  section  1,  the  concept  of  first-positivity  Is  Introduced,  and  some 
elementary  properties  developed.  In  section  2,  there  Is  a digression  on  some 
properties  of  relations  which  will  be  useful  In  the  sequel.  Section  3 studies 
first-nonzero  and  connectivity  properties  of  the  sequence  formed  by  the  powers  of 


a matrix;  at  this  stage,  there  Is  no  restriction  to  Metzler  matrices.  Section  4 
discusses  a property  of  matrix  exponentials,  which  are  used  In  expressing  the 
solution  to  a system  of  differential  equations.  In  section  5,  we  draw  the 
previous  results  together  In  application  to  the  first-positivity  and  connectivity 
properties  of  the  sequence  of  powers  of  a Metzler  matrix. 

The  results  to  which  the  earlier  sections  were  leading  are  contained  in 
section  6 and  section  7.  These  characterize  the  geneses  of  systems  of  difference 
equations  and  of  differential  equations,  respectively.  In  the  first,  the  time  of 
first-positivity  of  a specific  component  is  expressed  in  terms  of  the  connectivity 
properties  of  the  matrix,  the  specification  of  the  positive  components  of  the 
starting  values,  and  the  first-positivity  properties  of  the  forcing  function. 

In  the  second,  the  order  of  Increase  of  a specific  component  Is  expressed  in  terms 
of  the  same  factors. 

1.  First-positivity. 

In  the  following,  we  use  "D"  as  an  abbreviation  for  "Definition."  In  general, 
small  letters,  such  as  "x,"  stand  for  sequences  (functions  on  the  non-negative 
inters)  of  scalars  or  vectors;  a capital  letter,  such  as  "A,"  stands  for  a 
sequence  of  matrfiis;  or  "A(n)H  will  be  the  value  of  x or  A,  respectively, 

at  Integer  n;  "xa"  or  "Aag"  will  be  the  staiar  3jq"*oce  formed  by  considering  only 
the  component  a of  the  vector  or  the  element  (a, 3)  of  the  matrix  Sequtn-a. 

D.l.  For  a scalar  sequence  x,  define  v(x)  = min  ( n | x(n)  i 0 } . We  refer 
to  v (x)  as  the  non-zero  Index  of  x. 

Note  that  the  non-zero  Index  is  not  necessarily  defined,  since  x might  be 
Identically  zero.  This  problem  of  lack  of  definition  occurs  persistently  but  can 
*be  dealt  with,  as  will  be  noted  later.  * 

D.  2.  xv  $ x[v(x)] 

Thus  x Is  the  value  of  the  fir*;  clement  In  the  sequence  x which  Is  not  zero.  The 
symbv  l,  A"  Is  read,  "en*-^**  by  definition." 

I 3.  (x  Is  first  -positive)  ^(xv>0). 

Nov,  let  a be  a sequence  of  vectors.  Each  of  the  above  definitions  is  still 
allowed  tu  Sold,  but  must  be  applied  to  each  component.  That  Is,  for  each 
component  a of  the  vector,  xo  Is  a scalar  sequence.  Thenv(^)  is  defined  by  D.  1, 
and  similarly,  (xa)v  by  0.  2.  It  will  be  useful  to  Interpret  x to  be  a function 
from  a finite  domain  (the  domain  of  Its  components)  to  the  space  of  scalar 
sequences.  Let  F be  this  domain.  Then  v(x)  will  be  Interpreted  as  a function, 
the  value  for  a being  v(xa).  It  must  be  emphasized,  however,  that  \<x)  might  not 
be  defined  for  all  elements  of  F,  for  It  can  happen  that  the  sequence  xQ  Is 
Identically  zero  for  some  * and  therefore  v(x)  (o)  « v(^)  Is  not  defined  for 
that  value  of  a 

0.  4.  If  x is  a sequence  of  vectors,  the  function  v(x)  Is  defined  by  the 


/ 


relation, 
D.  5. 


v(x)(o) 
xv(a)  « 


‘ ^(xa). 
x«lv(xJ]  " x„ 


so  that  the  function  x\>  has  the  same  domain  of  definition  as  v(x). 

The  analogue  of  D.3  for  vectors  Is  a vector  Inequality,  l.e.,  we  regard  a 
sequence  of  vectors  as  first-positive  If  each  component  Is.  However,  we  have  to 
account  for  the  possibility  that  the  domain  of  definition  of  x might  not  be  the 
entire  possible  domain  F.  We  need  a convention  for  the  meaning  of  an  expression 
like  f(a)  > 0,  or,  more  generally,  for  an  expression,  f(a)  >g(a),  where  one  or 
both  of  the  functions  f and  g may  be  undefined  at  a particular  point  a.  We  shall 
regard  the  inequality  as  holding  If  f Is  undefined  ataand  g Is  defined  there  or 

If  both  are  defined  and  f(a)  is  Indeed  bigger  than  g(a),  and  not  otherwise. 

Convention  1.  If  the  functions  f and  g are  both  defined  at  a,  then  the 

expression,  f(a)  > g(a),  has  its  usual  meaning;  otherwise.  It  holds  If  and  only 
if  f Is  not  defined  at  a and  g Is. 

A similar  convention  will  hold  for  equality. 

Convention  2.  If  the  functions  f and  g are  both  define  at  a , then  the 
expression,  f(a)  * g(a),  has  its  usual  meaning;  otherwise.  It  holds  If  and  only 
If  neither  f nor  g are  defined  at  a. 

If  we  use  the  usual  symbol , 

0.6.  dom  f * {a  | f(a)  is  defined)  , 

then  Convention  1 says  that  f(a)  > g(a)  If  a e dom  g and  at  dom  f,  while 
Convention  2 Implies  that  f(a)  * g(a)  If  a t dom  f,  a i dam  g. 

Inequality  among  functions  has  the  usual  meaning  that  the  Inequality  holds  for 
all  values  of  the  argument,  with,  however,  the  conventions  above  observed.  From 
the  above  remarks.  It  Is  easy  to  note  that, 

f ■ g Implies  dom  f C dom  g. 

With  these  conventions,  the  definition  of  first-positivity,  D.  3,  remains 
valid  for  sequences  of  vectors.  It  means  that  each  component  Is  first-positive  If 
not  Identically  zero. 

The  following  simple  lemna  holds  for  adding  first-positive  sequences,  whether 
of  vectors  or  of  scalars. 

Irnmna  1.  If  x1  Is  first-positive  for  each  1,  then  E x1  Is  first-positive, 
and  v(e  x ) ■ min 


7-) 

1 

Proof: 


1 


v(x*). 


1 


mjn  v(x*),  P ■ { 1|  v(x*)  * v }.  If  n < v,  then  n<v  (x*) 
x*(n)  • 0 for  all  1,  by  0.1,  and  therefore  E x*(n)  • 0. 


Let  v 

for  all  1,  so  that  x'lnj  « o for  all  1,  by  0.1,  and  therefore  j: 

On  the  other  hand,  x^v)  • x1v  > 0 for  1 e P,  by  0.2  and  then  0.3,  vtx1)  > v for 
W P,  and  therefore  x1^)  • 0 for  1 4 p.  Hence,  E x1^)  > 0.  so  that  the  Lemma 
holds.  1 

Although  the  proof  has  been  stated  for  scalars.  It  holds,  with  suitable 
interpretation  for  vectors;  the  operator,  "min,"  In  the  statement  of  the  Lenina 
wist  be  taken  to  hold  component-wise. 


Also,  It  Is  useful  to  note  how  the  Conventions  are  used.  The  proof  as  given 
seems  to  require  that  v(x*)  Is  defined  for  all  1.  If,  however.  It  Is  defined 
only  for  some  1,  then  v Is  taken  as  the  minimum  over  all  1 for  which  It  Is 
defined.  That  Indeed  Is  the  Interpretation  Implied  by  Conventions  1 and  2.  Then 
the  argument  Is  valid  In  every  detail;  In  particular,  If  v(x*)  Is  not  defined  for 
some  1,  then  certainly,  by  Convention  1,  v(x^)  > v.  On  the  other  hand,  in  just 
that  case,  x*(n)  * 0 for  all  n,  and  therefore  certainly,  x^v)  = 0.  In  the  future, 
all  proofs  will  be  carried  on  as  though  all  functions  were  defined;  the  correction 
for  the  cases  of  lack  of  definition  can  easily  be  supplied  by  the  reader. 

The  same  definitions  will  be  needed  for  sequences  of  matrices;  however,  since 
a matrix  can  be  thought  of  as  a vector,  there  Is  no  need  for  additional  definitions. 
(Square)  matrices  can  be  thought  of  as  functions  on  a domain  of  the  form,  F x F, 
there  F Is  finite,  and  "x"  denotes  Cartesian  product.  In  this  case.  If  A Is  a 
sequence  of  matrices,  v(A)  Is  a function  of  two  variables,  representing  the  rows 
and  the  columns. 

A useful  concept  In  expressing  the  solution  of  systems  of  difference  equations 
Is  the  convolution  of  two  sequences,  a term  borrowed  from  probability  theory.  It 
is  the  same  as  the  expression  for  the  distribution  of  the  sum  of  two  Independent 
non-negative  random  variables. 

0.  7.  If  x and  y are  two  sequences,  then  the  sequence  x*y  Is  defined  by, 

(x*y)  (n)  - l x(j)  y(n-j). 

J-0 

First,  suppose  x and  y are  scalar  sequences.  From  D.l,  x(j)  ■ 0 If  j < v(x), 
and  y(n-j)  ■ 0 If  n-j  < v(y),  or,  equivalently,  j > n-v(y).  Hence, 
x(j)  y(n-J)  • 0 unless  v(x)  « .1  ■ n-  v(y). 

If  n < v(x)  + v(y),  then  v(x)  > n - v(y),  so  that  x(j)  y(n-j)  ■ 0 for  all  j,  and 
therefore  (x#y)  (n)  ■ 0.  If  n » v(x)  + v(y),  then  x(j)  y(n-j)  - 0 except  for 
J ■ v(x)  ■ n - v(y),  so  that  (x*y)  (n)  • x(v(x)j  y[v(y)l  ■ xvyv  f 0,  by  D.2  and 
D.l.  Hence,  for  scalar  sequences,  (x*y)  (n)  » 0 for  n < v(x)  + v(y),  f 0 for 
n ■ v(x)  ♦ v(y),  so  that,  by  0.1,  v(x*y)  * v(x)  + v(y),  and,  by  0.2,  (x*y)v  = 
xvyv.  If  It  Is  also  assumed  that  x and  y are  first-positive,  then,  by  0.3, 
xv  > 0,  yv  > 0,  and  therefore  (x*y)v  > 0,  so  that  x*y  Is  first-positive. 

\JSBLl-  If  x and  y are  scalar  sequences,  then  v(x*y)  » v(x)  ♦ v(y),  and 
(x*y)v  • xv  yv  . If  In  addition  x and  y are  first-positive,  then  so  Is  x*y. 

The  definition  of  a convolution  can  be  applied,  not  only  to  scalars  but  also 
to  vectors  and  pairs  consisting  of  matrices  and  vectors,  with  the  proper  Interpre- 
tation of  multiplication  In  the  definition.  Let  x and  y be  vector  sequences,  each 
vector  being  of  tht  same  number  of  components,  and  let  multiplication  be  Interpreted 
as  the  taking  of  an  Inner  product, 

*(J)  • i «j(J)  yj(n-j). 


so  that. 


(x*y)  (n)  • 2 E xfl(j)  yfl(n-J)  - EE  xa(J)  yft(n-j)  « E (x  *y«)  (n). 

j«o  0 8 0 ej-o  0 8 a 8 8 

Suppose  In  addition  that  x and  y are  first-positive.  Then,  for  each  0,  x0  and  yD 

M P 

are  first-positive,  and  therefore  Xg*y0  Is  first-positive,  by  Lemma  2,  and  x*y  Is 

first-positive  by  Lemma  1.  From  Lemmas  1 and  2, 

v(x*y)  - min  v(xg*y0  ) * min  [v(x)  (8)  + v(y)  (0)]  . (2) 

8 8 

He  Mill  apply  this  result  to  multiplication  of  a matrix  sequence  A by  a 
sequence  of  conforming  column  vectors,  x,  where  both  A and  x are  first-positive. 
Equation  (2)  can  be  applied,  with  each  row  of  A,  In  turn,  replacing  x,  and  x 
replacing  y.  The  sequence  A*x  Is  a sequence  of  column  vectors. 

Lemma  3.  Let  A be  a first-positive  sequence  of  matrices,  and  x a first- 
positive sequence  of  column  vectors  conforming  with  A.  Then  v(A*x)  (a)  - min 
[v(A)(a,8)  + v(y)  (8)J  and  A*x  Is  first-positive.  8 

2.  Relations,  their  Powers,  and  Chains. 

A relation  Is  simply  a set  of  ordered  pairs.  In  the  main  application  In  this 
paper,  the  relation,  C(A),  Is  defined  by  the  condition  that  a C(A)  8 Is  and  only  If 
Ao8  * °*  for  a given  matrix  A;  thus,  C(A)  Is  the  set  of  all  ordered  pairs  for 
which  this  condition  holds.  Typically,  then,  a relation,  R,  Is  some  subset  of 
F x F.  In  this  section,  however,  the  relation  R Is  arbitrary. 

A parti cular ^Interesting  relation  Is  the  Identity  relation,  E. 

D.8.  (o  E 0)  ■ (o  - 0) 

Following  Quine  (4),  p.  213,  the  relative  product  of  two  relations,  R & S, 

Is  defined  by, 

0.9.  (a  R|S  0)  > (for  some  y,  a R y and  y S 8). 

Like  any  other  form  of  multiplication,  the  relative  product  can  be  used  to  define 
the  powers  of  a relation  Inductively. 

0.10.  R°  - E.  R1*1  • Rn  | R. 

This  definition  can  be  given  a useful  alternative  form  by  Introducing  the 
concept  of  an  R-chaln.  A 

0.11.  (o  Is  an  R-chaln  of  length  n from  a to  0 ) * (a  is  a function  defined 
on  the  Integers  0,...,n,  o( 0)  > a,  o(n)  - 0,  and  o(1-1)  R o(1)  for  1 « 1,...,n) 

In  short,  an  R-chaln  Is  an  ordered  sequence  of  n+1  elements,  such  that  the  relation 
R holds  between  every  successive  pair. 

It  Is  Intuitively  obvious  and  can  easily  be  proved  by  Induction  that  the 
relation  R*1  holds  If  and  only  If  there  Is  an  R-chaln  of  length  n connecting  the  two 
elements. 

Lemma  4.  a Rn  0 If  and  only  If  there  exists  an  R-chaln  of  length  n from  a to 

B. 

For  a given  a and  0,  there  may  be  R-chalns  of  different  lengths  from  one  to 
the  other  (of  course.  It  Is  also  possible  that  there  are  no  R-chalns  of  any  length 
from  a to  0). 


D.  12.  (a  is  a shortest  R-chaln  from  a to  6)  ■ (for  some  n and  all  m n, 
o Is  an  R-chaln  of  length  n from  a to  8 and  there  Is  no  R-chaln  of  length  m fro* 
a to  8). 

Shortest  chains  have  a property  which  will  be  useful  later. 

Lemma  5.  A shortest  R-chain  is  a one-one  function. 

Proof:  We  seek  to  prove  that  If  o is  a shortest  R-chain  from  a to  8 and  if 
oO])  « o(i2),  then  i]  = i2>  Let  j.  = min  (1,,  i2),  h = max  (i1#  i?)  - U 
Clearly,  h * 0 by  definition;  we  seek  to  prove  that  h * 0,  or,  equivalently,  that 
h CO.  To  this  end,  we  construct  an  R-chain  of  length  n-h  from  o to  e.  where  o 
Is  an  R-chaln  of  length  n;  since  o is  a shortest  R-chaln,  It  follows  by  definition 


that  n-h  ■ n,  or  h * 0.  By  assumption, 

o(1_)  - o(1,  + h).  (3) 

Define  a function,  o'  on  the  Integers  0,...,n-h,  as  follows: 

o'(1)  * o(1),  0 • 1 < 1_,  (4a) 

* o(1+h),  1*1*  n-h.  (4b) 


If  i - 0,  then  from  (4b)  and  (3),  o'(0)  - o(h)  * o(0).  If  1 > 0,  then  o'(0)  * o(0) 
from  (4a),  so  that  o'(0)  * o(0)  In  either  case.  Also,  from  (4b),  o' (n-h)  = o (n). 

Since  o Is  an  R-chaln  of  length  n from  a to  3,  o(0)  * a and  o(n)  = b*  hence,  we 
have  shoMi  that  o'(0)  * a and  o' (n-h)  * B. 

To  show  that  o'  Is  an  R-chaln  of  length  n-h  from  a to  3,  It  remains,  by  D.ll, 
to  show  that  o'(l-l)  R o'(1),  1 * 1 * n-h.  If  1 < 1,  this  follows  Immediately  from 
(4a).  If  1 > 1,  then  1-1  * i;  by  (4b),  o'(l-l)  - o(1+h-l),  o’(1)  - o(1+h),  and, 
since  o Is  an  R-chaln,  o(1+h-l)  R o(i+h),  and  therefore  o'(l-l)  R o'(1). 

Finally,  let  1 - 1;  then  1 > 0.  In  this  case,  o' (1-1)  « o(1-1),  while 
o'(1)  ■ o(1+h)  - o(l  ♦ h)  - o(D  - o(1),  by  (4a),  (4b),  and  (3).  Since  o Is  an 
R-chaln,  o(1-l)  R o(1);  hence,  o' (1-1 ) R o' (1 ).  Therefore  all  the  conditions  of 
0.11  are  satisfied  for  o',  so  there  Is  an  R-chaln  of  length  n-h  from  a to  8,  and 
therefore  h * 0,  verifying  the  lemma. 

To  any  fixed  a,  there  Is  associated  the  set  (a  |oR  B)  ; we  Introduce  the 
notation  R B to  stand  for  this  set.  More  generally.  If  there  Is  a set  of  values  of 
6.  say  S,  then  R S,  will  stand  for  the  union  of  the  sets  RB,  for  6 c S. 

0.13.  R S • (a  | for  some  B,  a R e and  6 e S). 

It  Is  easy  to  see  that, 

R(SwT)«  (R  S)  v (R  T),  (5) 

where  Av8  Is  the  union  of  the  sets  A and  B. 

More  generally.  If  Sj  Is  an  index  set  of  sets,  with  1 varying  over  a set  1, 

«US, -y«$, . (6) 

where  V Is  the  union  of  the  sets  Sy 

Connectivity  of  a Matrix  and  Monzero  Index  for  the  Sequence  of  Powers  of  a 

Tht  following  connectivity  relation  can  be  associated  In  a natural  way  with 

| 1 

/ i 


« 


any  given  matrix  A. 

0.14.  (a  C(A)  6 ) - (AaB  f 0). 

The  connectivity  relation  for  a power  of  A will  also  be  of  Interest. 

The  sequence  of  powers  of  a matrix,  An,  n * 0,  where  A0  * I,  Is  a particular 
sequence,  and  a nonzero  index  function  can  be  associated  with  that  sequence.  As  a 
matter  of  notation,  we  must  distinguish  between  a particular  element  of  the 
sequence  and  the  name  of  the  sequence.  For  this  purpose,  we  borrow  the  functional 
abstractor  notation  from  mathematical  logic  (see,  e.g.  (4),  p.  226).  In  general, 
for  any  function  which  takes  on  the  value  f(x)  at  the  point  x,  we  will  mean  by 
*x  f(x)  the  name  of  the  function  which  takes  on  these  values.  In  this  paper,  we 
will  apply  the  notation  only  to  sequences,  where  the  variable  Is  n. 

D.15.  *n  f(n)  Is  the  function  which  takes  on  the  value  f(n)  when  the 
argument  takes  on  the  value  n. 

Thus,  xn  An  is  the  sequence  of  powers  of  A.  Associated  with  chls  sequence  is 
a nonzero  Index  function,  v(ar  An),  defined  over  ordered  pairs.  By  0.1,  0.4,  and 
0.14; 

v(An  An)  (o,  B)  * m1n{  n |a  C(An)e  >.  (7) 

Note  that, 

C(A°)  - C ( I ) - E.  (8) 

so  that  a C(A°)e  holds  If  and  only  If  o * e;  therefore, 

v(An  An)(a,e)  » 0 If  and  only  If  a » B.  (9) 

In  view  of  (7),  v (Ap  An)  Is  defined  for  a particular  pair  (o,b)  if  and  only  If 

aC(An)a  for  some  n.  In  set-theoretic  language,  define  the  relation,  Kv  (A),  by 

0.16.  Kv  (A)  - U C(An), 

• n«0 

where  U Is  the  union  of  the  relations  C(An);  remember  that  a relatlot  Is  a 
n*0 

particular  kind  of  set.  Then, 

Kw  (A)  ■ dom  v(An  An)  ■ dom  (Afl  An)v  . (10) 

He  will  explore  that  effects  on  the  nonzero  Index  and  related  concepts  for  a 

sequence  of  powers  of  a matrix  of  altering  the  matrix  by  adding  a constant  to  the 

dlogonal  elements.  As  a preliminary.  It  Is  noted  that  the  binomial  theorem  Is 

valid  for  pairs  ofnmatr1ces  which  commute  with  each  other. 

(A  ♦ B)"  - E (?)  Ar  Bn‘r  If  AB  ■ BA. 
r-0 

Let  B be  a scalar  multiple  of  the  Identity  matrix,  B ■ si  for  some  scalar  s.  Then 
si  commutes  with  every  matrix  A.  Note  that  (si )n_r  « $n-r  in-r  ■ *n“r  I. 

(A  ♦ sl)n  * E O s"-r  Ar  (11) 

fmQ  r 

0.17.  (A  » B mod  I)  ■ (A-B. « si  for  some  scalar  s) 

Note  that  the  relation  Is  symaetrlc.  Suppose  It  holds  for  two  matrices,  A and  B. 
From  (11),  n 

<B»U  ■ 1 (J)  •n*r  for  •**  scalar  s. 

r»0 


(12) 


Clearly,  If  (Bn^e  f 0,  then  It  must  be  that  (Ar)aB  f 0 for  some  r = n.  In 
particular,  let  n = v(An  Bn)(a,e).  Then  It  follows  that, 
v(An  An)(o,fl)  3<  v(An  Bn)(a,e). 

But  since  the  relation,  A * B mod  I,  Is  symmetric,  this  must  hold  with  A and  B. 
Interchanged;  also,  It  holds  for  any  a and  e. 
v(An  An)  « v(An  Bn). 

If  two  functions  are  equal,  they  have  the  same  domain  of  definition,  by  Convention 
2;  hence,  from  (10),  KV(A)  = Kv  (B).  Further,  if  we  set  n ■ v(An  An)(a,e)  = 
v(An  Bn)(a,e).  then,  by  (7),  (Br)  * 0 for  r < n;  from  (12),  (Bn)  = (An) 

By  0.5,  06  06 

(\,  An)v  = (An  Bn) v . 

Theorem  1.  If  A = B mod  I,  then  (a)  v(ar  An)  » v(a  Bn); 

(b)  (An  An)v  * (An  Bn)v  ; and  (c)  Kv  (A)  * Kv  (B).  " 

The  relation  C(A)  measured  what  might  be  termed  the  direct  connectivity  of  the 
matrix.  Two  elements  may  be  indirectly  connected  through  a chain  of  direct 
connections.  In  view  of  Lemma  4,  It  Is  natural  to  define  the  connectivity  index 
of  a matrix  (a  function,  not  a number)  as  smallest  power  of  C(A)  which  holds 
between  two  elements. 

0.18.  6(A)(a,8)  3 min  { n|  a(C(A)]  n 6 } 

Note  that,  since  [C(A)]°  - E,  6(A)(a,e)  * 0 If  and  only  If  a « 8. 

Also,  If  Ao6  f 0 and  a i s,  then  alCfMl’fi  while  not  a C(A)  °b,  so  that  6(A)  (a.e) 

* 1.  From  0.18,  the  domain  of  definition  of  the  connectivity  Index  is  precisely 

the  set  of  ordered  pairs  for  which  the  relation  (C(A)]  n holds  for  some  n.  Define, 

0.19.  K4  (A)  - U [C(A)J  n. 
n*0 

Then, 

K4  (A)  ■ dom  6 (A).  (13) 

It  Is  convenient  to  Introduce  the  notation, 

i(A)  (a.e  , n)  Is  the  set  of  C(A)-cha1ns  of  length  n from  a to  b.  (14) 

Then,  from 

a[C(A)]nB  If  and  only  If  z(A)  (a,  e,  n)  Is  non-empty.  (15) 

Me  Investigate  the  effect  on  the  connectivity  Index  of  a change  In  the 
diagonal  elements  of  the  matrix.  (Although  for  our  later  purposes,  only  a constant 
change  Is  relevant,  the  results  hold  for  any  change  In  the  diagonal  elements.) 
Suppose  therefore  A - B Is  a diagonal  matrix;  note  again  that  this  relation  Is 
symmetric.  Let  a be  a C(A)-cha1n  of  length  6 (A)(a,$).  By  0.18,  It  Is  a shortest 
C(A)-cha1n  from  a tv  and  therefore  a Is  a one-one  function  by  Lemma  5.  In 
particular,  o(1-l)  f o( 1),  1 * 1 * «(A)(o,e).  Since  A-B  Is  diagonal, 

**(1-1),a(1)  " Bo(1-1),  o(D*  06) 

Since  o Is  a C(A-cha1n,  o(1-l)  C(A)  o(1)  for  all  1,  or,  by  0.14,  A^^-j j 1 o, 
all  1.  By  (16),  t 0 lor  all  1,  so  that  o Is  a C(B)-cha1n.  In  the 


f 


notation  (14), 

£(A)  (a,  6,  5(A)  (a,e)  )C  Z(B)  (a,  8,  6 (A)  (<*,8)  ). 

If  6(A)  Is  indeed  defined  at  (o.s),  then  the  left-hand  set  Is  non-empty  and  there- 
fore so  Is  the  right-hand  set.  By  D. 18  and  (15), 

6(B)  (a, 6)  = 5(A)  (a, 8). 

This  statement  also  holds  if  6(A)  Is  undefined  at  (a,e)  by  Conventions  1 and  2.  By 
the  symmetry,  the  inequality  holds  with  A and  B interchanged,  so  that  6(A)  = 6(B). 
Since  the  two  functions  are  equal,  they  have  the  same  domain  of  definition.  From 
(13),  then,  we  have 

Theorem  2.  If  A - B is  a diagonal  matrix,  then  6(A)  * 6(B)  and  K6  (A)  * K6  (B). 
4.  Matrix  Exponentials 

As  Is  well  known,  solutions  to  systems  of  linear  differential  equations  with 
constant  coefficients  can  be  expressed  simply  in  terms  of  the  exponential  of  the 
matrix  of  coefficients. 

D.20.  e*  = z An/n! 

The  Infinite  seri8s°converges  absolutely  for  all  A,  so  that  eA  is  defined  and  the 

series  can  be  rearranged  at  will.  If  we  add  the  scalar  multiple  of  the  identity 

matrix  to  A,  the  value  of  the  exponential  can  be  expressed  with  the  aid  of  (11). 

eA+s1  = z (A+sI)n/n!  = z ( 1/n ! ) z (n!/l!j!)  s1  Aj 
n*0  n*0  1+j*n 


- Z Z (s1/H)(AJ/j!)  «(  Z sVl!)  ( Z AJ/j!) 
n-0  1+j»n  1*0  j*0 

= es  eA 

Lemma  6.  If  A * B mod  I,  then  eB  * p eA  for  some  positive  scalar  p. 

Corollary  1.  If  A * B mod  I,  then  eA  » 0 If  and  only  If  eB  ■ 0. 

Corollary  2.  If  * B mod  I,  then  C(eA)  * C(eB). 

5.  First-Positivity  and  Connectivity  of  Non-negative  and  Metzler  Matrices. 

First,  make  the  obvious  observation  that  If  A * 0,  then  aC(A)e  If  and  only  if 
A a > 0.  Recall, 


a 8 


(»"*')..  ■ M*">, 


„ '<*6  '«r  "VB 

If  A * 0,  then  An  * 0 for  all  n.  Hence,  the  right-hand  side  Is  a sum  of  non- 
negative terms  and  Is  positive  if  and  only  If  at  least  one  Is  positive.  Therefore, 

> 0,  or,  equivalently,  if  and 

07) 


oC(An+1)  8 If  and  only  If,  for  some  6,(An)aft  Aag 
only  If,  for  some  y,  oC(An)Y  and  >C(A)8  . In  the  notation  of  D.9, 

C(An+1)  • C(An)  | C(A). 

Lemma  7.  If  A « 0,  C(An)  - (C(A)]  n 

Proof:  For  n « 0,  we  know  that  C(A°)  « E • [C(A))°.  Suppose  the  Lemma  Is 
true  for  n.  Then,  from  (17), 

C(An+1)  • [C(A)]n  ( C(A)  - (C(A)ln+1, 


by  0.10. 

From  (7)  and  the  definition  of  6(A),  0.18,  Lemma  7 Immediately  Implies  that 


v( AR  An)  * 6(A)  when  A * 0;  the  two  functions  have  the  same  domain  of  definition, 

so  that  Kv  (A)  * K6  (A)  by  (10)  and  (13). 

Since  An  « 0,  it  follows  immediately  from  0.5  that  (xn  An)v  = 0;  but  by  0.4, 

It  must  be  that  (xfl  A0)''  (a,e)  t 0 for  all  (a,e)  In  the  domain  of  definition,  so 

that  (x„  An)v  > 0,  that  is,  the  matrix  sequence,  xn  An,  is  first-positive. 

Finally,  from  the  definition  of  an  exponential,  D.20,  e - 0 when  A * 0. 

Since  the  defininir.g  series  is  a sum  of  non-negative  terms, 

(eA)  a > 0 if  and  only  if  (An)  . > 0 for  some  n, 
ap  ap 

so  that, 

C(eA)  = U C(An)  = Kv  (A), 
n*0 

by  0.16. 

Lemma  8.  If  A » 0,  then  xfl  An  is  first-positive,  v(xn  An)  = 6(A),  eA  = 0,  and 
Kv  (A)  = K6  (A)  = C(eA). 

The  main  mathematical  result  of  the  paper  is  that  this  Lemma  holds  not  merely 
for  non-negative  but  for  all  Metzler  matrices.  We  recall  the  definition. 

0.21.  A is  Metzler  If  A . = 0 for  a 1 B. 

ap 

A simple  and  useful  relation  between  Metzler  and  non-negative  matrices  is  the 
foil  owl ng: 

A is  Metzler  If  and  only  If  there  exists  B * 0 such  that  A = B mod  I.  (18) 
Theorem  3.  If  A is  Metzler,  then  xn  An  is  first-positive,  v(xn  An)  = o(A), 
eA  « 0,  and  Kv  (A)  * K6  (A)  = C(eA). 

Proof:  Choose  B as  In  (18).  Then,  from  Theorem  1,  Lemma  8,  and  Theorem  2, 

v(*n  A")v*  v(An  Bn)v*  6(B)  = 6(A). 

From  Theorem  1 and  Lemma  8, 

(X„  AV  - (Xn  bV  > 0. 

From  Lemma  8,  eB  * 0,  and  therefore  from  Corollary  1,  eA  = 0.  Finally,  from 
Corollary  2,  Lemma  8,  and  Theorem  1,  C(eA)  - C(eB)  = Kv  (B)  * Kv  (A),  while  from 
Lenina  8 and  Theorem  2,  Kv  (B)  « K6  (B)  * Kv  (A). 

The  Importance  of  this  theorem  Is  that  the  qualitative  behavior  of  the  powers 
and  the  exponential  of  a Metzler  matrix  can  be  Inferred  solely  from  its  connectivity 
properties.  These  depend  only  on  the  location  of  the  off-diagonal  zeroes  and  are 
Independent  both  of  the  diagonal  elements  and  of  the  magnitudes  of  the  non-zero 
off-diagonal  elements.  Thus,  If  we  raise  a Metzler  matrix  to  successively  higher 
powers  we  know  that  In  each  place  in  the  matrix,  the  first  non-zero  element  (if 
any)  will  be  positive  and  the  power  for  which  the  non-zero  entry  occurs  Is  equal  to 
the  length  of  the  shortest  chain  from  the  row  element  to  the  column  element  through 
non-zero  entries. 

A side  consequence  of  the  analysis  Is  a pair  of  what  are  apparently  new 
necessary  and  sufficient  conditions  for  a matrix  to  have  the  Metzler  property. 

Theorem  4.  Each  of  the  following  conditions  Is  necessary  and  sufficient  that 
A be  a Metzler  matrix:  (a)  xn  An  Is  first-positive;  (b)  eAt  ■ 0 for  all  t > 0. 


Proof:  (a)  Necessity  has  already  been  shown  In  Theorem  3.  Suppose,  then,  that 

.n 


the  sequence  xR  A"  Is  first-positive.  The  pairs  (0,0)  can  be  classified  according 

as  v( A An)  (0,0)  Is  0,  1,  or  greater  than  1.  In  the  first  case,  as  remarked  In 
n 1 

(9),  a = 0.  In  the  second,  we  must  have  (A  ) f 0,  by  definition,  and  therefore 

A > 0,  since  A An  Is  first-positive.  In  the  third  case,  (A^ ) „ * 0 by 

op  n ^ CL  p 

definition  of  the  non-zero  index.  Hence,  Ifaf 
Metzler  matrix. 

(b)  If  A Is  Metzler  and  t > 0,  a scalar,  then  At  Is  also  Metzler,  and  eAt  = 0 


> a6 
6,  A * 0,  so  that  A Is  a 


At  ^ 

by  Theorem  3.  Conversely,  suppose  that  e * 0 for  all  t > 0. 
when  t * 0,  that. 


Note  that  eAt  = I. 


d(eAt) 


^t 
so  that, 
d 


A e 


At 


'P 

|t=( 


A, 


and  that,  by  definition, 
d W 


dt  U 


11m 

t-0 


ii. 


From  the  hypothesis  and  (18),  eAt  - I is  Metzler  for  t > 0,  and  therefore 
(eAt  - I)/t  Is  a Metzler  matrix  for  t > 0.  Hence,  A Is  a limit  of  Metzler  matrices; 
since  the  set  of  Metzler  matrices  Is  clearly  closed,  from  the  definition,  A must  be 
a Metzler  matrix. 

6.  Genesis  of  a System  of  Difference  Equations  with  Metzler  Matrix  as  Jacobian 
We  consider  the  system  of  difference  equations, 

x(n+l)  « Ax(n)  + b(n),  (19) 

where  we  assume  that  A Is  a Metzler  matrix,  b a first-positive  sequence,  and 
x(0)  * 0.  To  express  the  solution  compactly,  define  a vector  sequence,  c,  by, 
c(n)  ■ x(0)  for  n * 0, 

* b(n-l)  for  n > 0.  (20) 

Then  (19)  can  be  written, 

x(n+l)  ■ Ax(n)  + c(n+l), 
x(0)  * c(0). 

By  Induction,  It  Is  easy  to  verify  that, 

x(n)  * ? A^  c(n-j), 

J-0 

or.  Ip  the  notation  Introduced  In  0.7, 

x-(AnAn)*c.  (21) 

Since  b Is  first-positive.  It  Is  easy  to  see  from  (20)  that  c Is  first-positive  (If 
xg  (0)  > 0,  then  cg  (0)  > 0,  so  that  cg  Is  certainly  first-positive;  If  xg  (0)  * 0, 
then  the  first  non-zero  element  In  the  sequence  cg  Is  the  first  non-zero  element 
of  b.  with  non-zero  Index  Increased  by  1,  and  must  be  positive  since  b Is  first- 


positive.)  If  A Is  Metzler, 


A A 
n 


Is  first-positive  by  Theorem  3.  Hence,  by 


(22) 


Lenina  3,  x Is  first-positive,  and, 

v(x)(o)  - min  l v(ar  An)(a,e)  + v(c)(b) J . 

Define,  for  any  vector  x, 

D.22.  P(x)  $ {«  |x  > 0 }. 

In  (22),  for  each  s,  either  6eP[x(0)j  or  8/  P[x(0)]  . The  minimum  can  be  taken 
separately  over  the  two  sub-sets  and  then  the  minimum  of  the  two  taken.  Note  that 
If  6 eP[x(0)],  then,  from  ( 20),  v(c)(e)  = 0,  while  If  B i P[x(0)j,  then  v(c)  (e)  = 
v(b)(e)  + 1.  Further,  from  Theorem  3,  v(xfl  An)(o,0)  ■ 6(A)(a,0).  Substitution 
Into  (22)  yields. 

Theorem  5.  If  A Is  a Metzler  matrix,  b a first-positive  sequence  of  vectors, 
x(n+l)  » Ax(n)  + b(n),  and  x(0)  * 0,  then  x Is  a first-positive  sequence,  and 
v(x)(a) 


« min  { min  6(A)(a,e),  1+  min  [«( A)(a,e)  + v(b)(e)]}  . 
6 cPlx(O)]  e *P[x(0)] 


Note  that  Theorem  5 Implies  that  each  component  is  positive  before  it  can 
become  negative.  Further,  a given  component  can  be  positive  In  two  different  ways. 
One  Is  ultimately  due  to  a positive  Initial  component  a which  Is  linked  to  the 
given  component  a directly  or  Indirectly.  The  other  Is  through  the  emergence  of 
a positive  element  In  one  component  of  the  forcing  term  b(n),  which  Is  then  linked 
to  the  given  component,  o.  The  shortest  of  all  these  routes  determines  the  length 
of  time  before  the  positive  effect  appears. 

7.  Genesis  of  a System  of  Differential  Equations  with  Metzler  Matrix  as  Jacobian 

As  a preliminary,  we  note.  In  the  notation  Introduced  in  D. 13  and  0.22, 

Lewwa  9.  If  A - 0 and  x * 0,  then  P(Ax)  - C(A)  P(x). 

Proof:  (A^)  ■ | A^  xg  . Since  all  terms  are  non-negative  by  assumption, 
(Ax)a  > 0 If  and  only  If,  for  some  6,  A^  > 0 and  xg  > 0;  but  this  holds  If  and 
only  If,  for  some  e,  aC(A)e  and  s eP(x). 

Now  consider  the  system  of  differential  equations, 

x - A ♦ b(t),  (23) 

where  A Is  a Metzler  matrix,  b(t)  * 0 for  all  t,  and  x(0)  ■ 0.  This  clearly  has 
the  solution, 

x(t)  - eAt  x(0)  ♦ fl  eA(t'u)  b(u)  du, 
o 


■ y(t)  ♦ z(t) , 

where, 

y(t)  ■ eAt  x(0), 

*(t)  • / «A(t'a)  b(u ) 
o 

By  Theorem  4.  eAt  • 0 
t ^ o.  Also,  eA't’u^  ■ 0 

eA(t"u)  b(u)  t o,  0 * 


(24) 


(25) 

du.  (26) 

for  all  t > 0,  so  that,  from  (25),  y(t)  * 0 for  all 
for  t > u;  since  b(u)  ■ 0 for  all  u,  by  assumption, 

u < t. 


(27) 


and  therefore  z(t)  0 from  (26).  Combining  these  statements  we  see  that, 

x(t)  3 0,  and  P(x(t)J  - P[y(t))v/  P[z(t)j  , all  t>0.  (28) 

From  lemma  9,  P{y(t) 1 3 C(eAt)  P(x(0)J,  all  t>  0. 

From  Theorem  3,  C(eAt)  ■ K 6 (At)  for  t > 0;  but  since  obviously  C(At)  * C (A) 
for  any  matrix  A and  any  scalar  t i 0,  It  follows  from  the  definition  of  K 6 (A), 

D. 19,  that  K 6 (At)  * K 6 (A)  for  all  t>  0. 

P[y(t) ] 3 K 6 (A)  P [x (0) ] for  t > 0.  (29) 

In  particular,  P[y(t)]  Is  independent  of  t for  t > 0.  Hence,  as  far  as  the  effects 
of  Initial  conditions  go,  any  component  that  Is  going  to  be  positive  eventually 
Is  positive  Immediately.  However,  as  will  be  seen  below,  the  delay  effects  of 
the  connectivity  of  the  matrix  affect  the  solution  but  in  a different  way. 

If,  for  some  component  a,(eA^t-u)  b(u)  ) * 0 for  all  u,  0 * u < t,  then 

obviously,  from  (26),  (z(t))a  3 0.  On  the  other  hand,  if  b is  assumed  continuous, 
then  if 

(eA't'u)  b(u))a  > 0 for  some  u,  0 * u<t  , 

It  is  positive  in  some  interval  and  hence  from  (27)  and  (26),  (z(t))Q  > 0.  In 
symbols, 

P[z(t)J  3 U P[eA(t"u)  b(u)]  3 U C(eA*t-u))  P[b(u)] 

03  u < t 0»  u < t 

3 < K6  (A)  P[b(u)] 

0-  u < t 

- K*  (A)  „ (J  P [b(u ) ] , for  t>0.  (30) 

03  u < t 

The  steps  are  the  same  as  those  leading  to  (29),  together  with  a final  step  which 
uses  (6). 

The  result  can  be  made  still  more  transparent  with  the  aid  of  a nonzero  Index 
for  functions  of  a continuous  variable. 

0.23.  c(b)  (8)  3 Inf  {t  | t 3 0.  bp  (t)  i 0}. 

If  u <c(b)  (b),  then  bfl  (u)  3 0,  and  therefore  B i P(b(u)].  Therefore,  If 
t • c(b)(B),  B t P(b(u) ] for  all  u,  0 3 u < t,  and  therefore. 

If  t 2 c (b)(8),  B t <U  P(b(u)J. 

Oh<t 

Suppose  now  t > c (b)(8).  Then,  by  0.23,  there  exists  u, 

c(b)  (6)  3 u < t,  8 eP[b(u)), 
and  therefore, 


c .U  P(b(u)) 
03  u < 


t 


Hence, 

U P[b(u) 1 If  and  only  If 
8 e 0*u  < t 

0 ^ ((b)(8)  < t.  (31) 

For  fixed  b,  c(b)  1$  a function  over  a finite  set.  For  any  function,  f, 
f’^(y)  swans  the  set  (x|  f(x)  ■ y);  for  any  set  S In  the  range  of  f,  f 1 (S)  means 


the  set  (x  |f(x)  e S).  In  this  notation, 

0 ■ c(b)(e)  < t If  and  only  If  BeU(b)]"1  (<0,t)), 
where  <0,t)  Is  the  Interval  closed  on  the  left  and  open  on  the  right. 

Then, 

< U P b(u)  Mc(b)]'1  (<0,  t)). 

0»  u < t 

In  combination  with  (30),  (29),  and  (28),  we  can  state. 

Theorem  6.*  Suppose  A Is  a Metzler  matrix,  b continuous  and  non-negative, 
x(t)  » Ax(t)  + b(t)  for  all  t * 0,  and  x(0)  * 0.  Then,  x(t)  * 0,  all  t = 0,  and, 

P(x(t)]  * K 6 (A)  (P[x(0)  o [C(b)j  'Vo.t)  )J  for  t > 0. 

Notice  again  that  the  signs  of  all  the  components  are  completely  determined 
by  the  sign  patterns  of  the  matrix  of  the  initial  conditions,  A and  of  the  forcing 
function  b.  At  any  time  t,  we  find  all  components  which  are  either  positive  at 
time  0 or  have  been  nonzero  at  some  time  point  before  t;  then  take  all  components 

linked  to  them  directly  or  Indirectly  by  chains  of  nonzero  entries  In  the  matrix 

A.  This  set  Is  precisely  the  set  of  positive  components  at  time. 

This  theorem  actually  relates  to  more  than  ttegenesis  of  the  dynamic  system. 
The  next  result  will  study  behavior  at  the  starting-point,  specifically,  the 
qualitative  behavior  of  the  successive  time  derivatives  of  the  different  components 
of  x(t)  at  the  point  t * 0. 

First,  the  non-negativity  of  a function  and  the  first-positivity  of  the 
sequence  of  its  derivatives  are  related. 

Lewwa  10.  An  Infinitely  differentiable  vector  function  f Is  non-negative  for 
t • 0 If  and  only  If,  for  each  t * 0,  the  sequence  xn  f^(t)  Is  first-positive. 

Proof:  Suppose  f(t)  * 0,  all  t * 0,  but  for  some  tQ  * 0,  the  sequence 
An  f 'V  is  not  f1rst"P°sit1ve*  Then  there  exist  s and  n so  that, 

^ (t0)  » 0 for  r < n,  f*n)  (tQ)  < 0.  But  If  n * 0,  then  f0  (tQ)  < 0, 
contrary  to  hypothesis;  If  n > 0,  then  f (t)  < 0 In  some  right-hand  neighborhood 
of  tQ,  again  contrary  to  hypothesis. 

Conversely,  suppose  *n  f*n*(t)  Is  first-positive  for  all  t.  Then  In 
particular.  It  Is  Impossible  that  fg  (t)  < 0 for  any  t and  6,  for  then  f^(t)  < 0, 
In  which  case  the  sequence  f^(t)  would  not  be  first-positive. 

Differentiate  the  system  of  differential  equations  (23)  n times,  and  then 
set  t ■ 0. 

x(n+1)(o)  . Ax^(0)  + b^(0).  (32) 

By  Lemma  10,  the  sequence  b(n)(0)  Is  first-positive;  the  matrix  A Is  Metzler 
by  assumption;  and  x^(0)  * x(O)  ■ 0,  by  assumption.  Hence,  (32)  constitutes  a 
system  of  difference  equations  which  satisfies  all  the  hypotheses  of  Theorem  5. 

Theorem  7.  Under  the  hypotheses  of  Theorem  6,  the  sequence  xnx^(0)  is 
first  positive,  and 

* The  conclusion  that  x(tj  * 0 already  appeared  In  (1),  Theorem  *,  p.  14. 


v(*n  x(n)(0)  ) 

-min{  min  6 (A)(a,e),  1 * min  6[(A)(a,e)  ♦ v(x  b(n,(0)l) 

B c P [x (0) J 8 i P[x(0)J  n 

Thus,  a positive  Initial  component  causes  every  component  Indirectly 
connected  to  it  to  become  positive  in  the  right-hand  neighborhood  of  the  origin, 
but  the  order  of  growth  (linear,  quadratic,  or  whatever)  depends  on  the  length 
of  the  connecting  chain  through  the  matrix.  Similarly,  a forcing  term  will  cause 
an  order  of  growth  In  a component  of  x which  Is  greater  by  1 than  the  sum  of  the 
order  of  growth  of  the  forcing  term  at  zero  and  the  length  of  the  shortest  chain 
to  the  x-component.  These  remarks  are  only  valid  for  the  first  effect  on  the 
given  component. 

To  Illustrate,  for  an  x-component  which  Is  Initially  zero,  the  growth  is 
linear  if  either  there  Is  a chain  of  length  1 to  a positive  x-component  or  the 
forcing  term  for  the  given  component  Is  positive.  The  growth  Is  quadratic  if 
neither  of  these  conditions  hold  and  If  one  of  the  following  three  conditions  is 
valid:  (1)  there  Is  a chain  of  length  2 to  a positive  initial  component;  (2) 
there  Is  a chain  of  length  1 to  a component  whose  forcing  term  Is  Increasing 
linearly  from  zero;  (3)  the  forcing  term  for  the  given  component  Is  increasing 
qjadratlcally. 
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This  paper  studies  the  behavior  In  the  neighborhood  of  >the  starting  point  of  dynamic 
systems  (solutions  of  difference  or  differential  equations)  whose  Jacobians  are  Metzler 
matrices.  (A  Metzler  matrix  is  one  whose  off-diagonal  elements  are  non-negative.) 

A new  concept,  that  of  first-positivity  of  a sequence,  is  introduced;  a sequence  is 
first-positive  if  its  first  non-zero  element  is  positive.  First-positivity  holds 
for  a sequence  of  matrices  or  vectors  if  it  holds  for  each  component.  It  is  shown 
that  the  sequence  of  powers  of  a Metzler  matrix  is  first-positive;  also,  for  each 
position  in  the  matrix  (defined  by  row  and  column),  the  number  of  steps  to  the  first 
non-zero  entry  in  the  sequence  is  equal  to  the  minimum  length  of  a chain  .from  the 
row  to  the  column  through  non-zero  entries  in  the  Metzler  matrix.  From  this,  it  is 
possible  to  express  (a)  the  time  of  flrst-posltlvlty  of  a specific  component  of  the 
solution  to  a system  of  difference  equations  and  (b)  the  order  of  Increase  of  a 
specific  component  of  the  solution  to  a system  of  differential  equations  in  terms  of 
the  connectivity  properties  of  the  governing  matrix,  the  specification  of  the  positive 
components  of  the  starting  values,  and  the  flrst-posltlvlty  properties  of  the  forcing 
function. 
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